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Abstract 

If the baryon asymmetry of the universe was created at the GUT scale, sphalerons together 
with exotic sources of (5— L)-violation could have erased it, unless the latter satisfy stringent 
bounds. We elaborate on how the small Yukawa coupling of the electron drastically weakens 
previous estimates of these bounds. 



1 Introduction 

The first proposals to explain the baryon asymmetry of the universe suggested that baryoge- 
nesis occurred at the very high temperatures characteristic of the scale of grand unification. 
While there has been great interest in generating baryons at the electroweak scale, the 
feasibility of this approach remains uncertain, whereas baryon-production in grand unified 
theories can be calculated fairly reliably in the framework of a given GUT. However there are 
many potential opportunities for a baryon asymmetry produced at 10^^ GeV to be destroyed 
later on, which leads to interesting constraints on the coefficients of operators that violate 
baryon or lepton number. 

A baryon asymmetry generated at the GUT scale will survive only if baryon number 
violating interactions are subsequently out of thermal equilibrium. In most SU{5) grand 
unification schemes, any baryon asymmetry generated will be accompanied by a lepton 
asymmetry such that B — L = 0. Non-perturbative electroweak interactions via sphalerons 
also violate baryon and lepton (actually B + L) number at high temperatures |I| while 
perserving B — L. As such, the GUT-scale produced baryon asymmetry will be erased if 
sphaleron interactions are in thermal equilibrium, as they are expected to be at temperatures 
between ~ 10^ and 10^^ GeV. By the same token, any GUT-scale baryon asymmetry with a 
net B — L asymmetry will be preserved. 

If in addition to the sphaleron interactions there was another baryon or lepton number 
violating interaction in a linear combination other than B-\-L and in thermal equilibrium, all 
net baryon and lepton numbers would be erased. Such was the case encountered by Fukugita 
and Yanagida (FY) |^ when considering a simple extension to the standard electroweak 
model to allow for Majorana neutrino masses. In order to maintain a net baryon asymmetry, 
this source of lepton number violation is constrained to be out of equilibrium at the times 
when sphalerons are in thermal equilibrium. This constraint was then translated into a 
bound on the neutrino mass. 

FY derived the conservative bound m^>j50 keV, assuming only that sphalerons are in 
equilibrium at T ~ 10^ GeV. It was subsequently argued that this constratint could be 
made substantially stronger because sphalerons are in equilibrium at significantly higher 



temperatures as well [Q. These arguments were also extended to other sources of baryon 
and lepton number violation p, |^ . 

In a previous paper [^ we pointed out that the constraints are dramatically weakened 
when one realizes that right-handed electrons can protect the baryon asymmetry in cases 
where it would otherwise have been destroyed, so long as the interactions of the operators 
in question fall out of equilibrium before the interactions of cr with ei and Higgs bosons 
come into equilibrium. The idea that the weakness of e/j-CL-Higgs interactions might have 
something to do with preserving the baryon asymmetry from sphalerons and other B- or 
L- violating effects was originally considered by the authors of ref. ^j (CDEO), for the initial 
condition that B — L = 0. Their goal was to determine whether cr ^-> ei equilibrium might 
be established so late that sphaleron interactions would not have enough time to nullify the 
baryon asymmetry. 

The argument for the importance of cr decoupling is intuitively very simple and was 
presented in a rudimentary form in our previous work [0. Since we feel that the idea 
has far-reaching consequences, it seems worthwhile to rederive the earlier results in a more 
accurate and rigorous fashion. Moreover there is a variety of special cases for the initial 
conditions or the depletion at intermediate scales of the B and L asymmetries, each with its 
own subtleties, the discussion of which was beyond the scope of the first paper. Filling these 
gaps is the goal of the present work. 

In what follows we will develop the full network of Boltzmann equations necessary to 
solve for the baryon asymmetry through the thermalization of cr. In section 3, we derive 
the rates for decays and inverse decays as well as scatterings to determine the cr decoupling 
temperature. We then consider the case with B — L = as did CDEO, in section 4. 
Although their result indicated that the baryon asymmetry was not preserved in this case, 
the final baryon asymmetry is exponentially sensitive to changes in reaction rates, and we 
have undertaken a careful reanalysis of some effects which were treated too roughly by 
CDEO. In particular, we have included the effects of scatterings. Qualitatively, our results 
confirm those earier conclusions. In section 5, we generalize to the case where B — L ^ and 
describe several mechanisms in which a baryon asymmetry can be preserved from sphaleron 
effects. Our conclusions will be given in section 6. 



2 The network of Boltzmann equations 

In this paper we explore the consequences of the fact that the right-handed chiral elec- 
trons cr remain out of equilibrium to suprisingly low temperatures in the early universe. 
This phenomenon arises simply due to the smallness of the electron mass or its Yukawa 
coupling, he — 2.94 x 10^^. At temperatures above the electroweak phase transition elec- 
trons are massless, or more accurately, different chiral states have only chirality-conserving 
temperature-dependent masses, so that the only processes that change lepton chirality are 
interactions with Higgs bosons.Q The lepton chirality-changing rate is therefore proportional 
to the small number h^. The smallness of the cr-cl equilibration temperature T^, is im- 
portant because the equilibrium conditions in the primieval plasma are different above and 
below T*. This has far-reaching consequences for the regeneration of the baryon asymmetry 
and leads to a considerable relaxation of bounds imposed on various exotic interactions to 
preserve the primordial baryon asymmetry. 

It was noted in CDEO that in the particular case of a universe with initial B — L = 
the final baryon asymmetry is exponentially sensitive to the value of the rate of the chirality 
changing interactions. In such a case special care must be taken in computing the thermally 
averaged rates and also in deriving the evolution equations for the asymmetries from the 
Boltzmann equations for the particle densities. Indeed the evolution equations found in the 
literature are often somewhat inaccurate in their treatment of the relevant collision terms in 
the Boltzmann equations. 

In this section we will carefully derive evolution equations for the particle asymmetries 
in the early universe during the time when cr ^->- cl transitions come into equilibrium. In 
general these equations couple all particle species present in the heat bath and involve a 
very large number of reaction rates. Fortunately a natural hierarchy exists between different 
reactions; various interaction processes can be divided into 'fast' or 'slow' depending on 
how their rates compare to the Hubble expansion rate of the universe H at the timescale 



^ For the quarks, QCD instanton interactions also violate chirality, but it is easy to show that these give 
no new conditions on the chemical potentials; they are consistent with the constraints associated with Higgs 
boson interactions. 



(Inverse) Decays 


H* ^ enL/ 


Gauge scattering 


CrH ^ LeA 


CrA ^ LeH* 


crLJ" ^ AH* 


Fermion scattering 


CrLj ^ LeJR 


CRLf^ ^ LejR 


CRLe" ^ LfjR 


Electron lepton 
number violation 


eiLi ^^ h^H 



Table 1: Shown are all tree level interactions that change the chirality of electrons and also a 
particular type of electron lepton number violation. H is the Higgs doublet, with components 
{h~^, hY . Li refers to either of the members in a left-handed doublet; / represents all fermions 
but a is restricted to the leptons. A can be either the B or the W gauge field. 



of interest. A process is called fast only if its rate is much bigger than H, otherwise it 
is considered to be slow. Fast processes need not be treated in detail, but instead, their 
effect can be accounted for by an appropriate set of equilibrium conditions between chemical 
potentials and certain conservation laws that restrict the explicit equations governing the 
slow processes. 

Examples of fast processes are the various electroweak interactions connecting fermions 
in the weak isospin doublets, W~ <-> i7eG~,ud, etc., and the interactions mediated by the 
electroweak anomaly, i.e. the sphaleron interactions [jl[. Naturally, all interactions that 
change the chirality of electrons are taken to be slow, but there can also be others, due 
to new physics. In Table 1 we have collected all slow processes that are relevant to our 
considerations in this section. For illustrative purposies we also show an electron lepton 
number violating interaction, which we will study in detail in section 5 below. 

Evolution equations for the asymmetries in the cr and ei species can be written down in 
a quite general manner, but to make the equations close typically requires some additional 
constraint that is specific to a given system. The problem shows up, for example, in the 
appearance of the Higgs chemical potential /iq in the evolution equations; the expression 



for /io in terms of fiea and /ig^, if it can be given, depends on the particular model. It is 
also clear that if any particle number other than ei or cr numbers was softly violated by 
some new exotic slow interaction, then additional equations governing the evolution of the 
corresponding asymmetry would be needed. Here we will restrict ourselves to the case where 
at most ei violation (in addition to cr-cl equilibration) is tracked in detail, which includes 
the special situations that we study in sections 4 and 5. Our derivation however will address 
all the subtleties needed to derive evolution equations for any other system of interest. 

2.1 Equilibrium conditions 

Let us first consider some of the equilibrium conditions established by the fast processes. 
Because sphaleron interactions fall out of equilibrium very quickly below the electoweak 
phase transition temperature ^ and because the sphaleron reprocessing of the changing eji 
asymmetry is crucial for our considerations in this paper, it is sufficient for us to consider the 
evolution equations only at temperatures above the phase transition. In the high temperature 
phase, the electroweak interactions mediated by the charged gauge bosons and scalars force 
the chemical potentials of quarks and leptons to obey the equilibrium conditions [^ 



;i) 



Here the labels u and d refer to all quarks with weak isospin +1/2 and —1/2 respectively, 

and i refers to the muon and tau lepton; the chirality equilibrium equation is not assumed 

to hold for electrons. fiQ is the chemical potential of the Higgs field. We have assumed 

that the universe was initially isospin-neutral,^ implying that the W chemical potential, 

^In contrast to the charge neutrahty of the Universe, which is proved to a very high degree by observations, 
there is no requirement for the universe initially to have had zero net isospin Q3; any primordial isospin 
would have vanished immediately after the phase transition leaving no observable trace today. However, it 
is improbable that such an asymmetry could have been dynamically generated because Q3 corresponds to a 
gauge symmetry. Moreover, any initial value for it would have been washed away by inflation. We therefore 
adopt the value Q3 = as the most natural choice, albeit not an absolutely necessary. 



/^di = /^«L 




/^ei = f^ue 


^J'UR — ^J'UL = /^o 




/i^i = fJ'ue 


fJ-dR - l^dL = ~/^0 


f^iR 


- /^£l = -/^O- 



which would otherwise appear in eq. ([1|) above, is zero at high temperatures 0]. Note 
also that due to the flavour mixing in the quark sector all quark flavours have the same 
chemical potential. Equilibrium sphaleron interactions impose a further condition between 
the chemical potentials of quarks and leptons 

9^UL + /^ei + Yl f^^L = 0> (2) 

e=ii,T 

and the charge neutrality of the universe leads to 

Q = QfJ'UL + 13/io - ficL - Ate_R - 2 Yl f^<^L = 0- (3) 

i=tl,T 

Equations d^-^) contain six unknown parameters. In addition to (@-^ we will always have 
four additional conditions on the chemical potentials so that the system has a unique solution. 
Assuming that there were no new interactions, then well above T* we would have the following 
conservation laws:0 



l^ 


-L, -- 


- ^e,p 


5- 


-L, - 


- C 


Lb- 

3 


-Lr -- 


- a,p 




^eR = 


- c 



(4) 

where Lf = Lf^ + L/^ + L^j and the subscripts p in the constants C indicate that they are 
the primordial values. The connection between the asymmetries and the chemical potentials 
is 

i/-^^-j^f + o((^,/rf). (5) 

where c is 1 for fermions and 2 for bosons and s = '^g^{T)T^ is the entropy density of the 
universe. 



■^We warn the reader that the symbol Le can mean the electron lepton doublet field, the asymmetry in 
electron lepton number, or the U(l) symmetry for electron lepton number, depending on the context. A 
similar warning applies for L^ and L^. 



If any of the conservation laws (^ are violated by exotic fast interactions, then they must 
be replaced by new equilibrium conditions. Indeed, the situation well below T* provides an 
example of how this works. There the e^ — cr transitions violate cr number strongly, so 
that the last equation in (H) is replaced by the equilibrium condition /ig^ — /ie^ = — /Uo- On 
the other hand, if some of the conservation laws are softly violated, so that around T^, the 
rates of relevant interactions are comparable to the expansion rate, then explicit equations 
must be written for the asymmetries affected by these interactions. In this section we will 
concentrate on the cases where Lg^ is violated by itself or in combination with soft electron 
lepton number violation. 

2.2 The evolution equations 

To derive the evolution equations for the e^j and e^ asymmetries we first write down the 
Boltzmann equations for the difference between the particle and antiparticle number den- 
sities. The collision terms in these equations correspond to the interactions shown in the 
table 1 and depicted in figure 1. It is reasonable to assume that each particle is in kinetic 
equilibrium at all times of interest; this is true in particular for the right-handed electrons 
due to their nonvanishing hypercharge. We then find 

(dt - 3H)in,^ - n,^c) = 2 J l[dn,i2nY S\p, - P2 - P3)\Mih* ^ e^'^en)]' 

1=1 

X {/,,(! - /,,0(1 - fen) - /e.^/ejl + M} " CC. 

+2ff[dn, (27r)V(pi +P2 -P2 -P3) E {iMihjR ^ CL^en)]' 
j=i / 

+ ---}-c.c. (6) 



+y'ndn, C2T,Yb\vx^V2-V2-Vz) E {Si\M{h-h- ^ hCL) 



i=l i=e,fj.,T 



+ ■■■}- c.c. + ^{dt - 3H){nB - riBc) (7) 

where H is the Hubble expansion rate and we defined the usual phase space volume element 
as dn = d'^p/{27T)^2E. J'l and fn are the chiral quarks or leptons. The symmetry factor Si in 
(6) is introduced to avoid double-counting in case of identical initial or final state particles: 
Se = \ and Sfj,^r = |- According to our assumption of prevailing kinetic equilibrium, the 
distribution functions are the usual Fermi-Dirac or Bose-Einstein functions 

/,(p,/.,) = (e^(^-^»)±l)-\ (8) 

where (3 = 1/T and T is the common temperature. The factors of 2 in front of the integrals 
in (^) come from summing over the equal contributions from both fields in each left handed 
doublet. This is evident because the couplings and therefore the temperature-dependent 
masses of the doublet components are identical. Similarly, the half in front of the first term 
in (0) arises because CL-number is changed only in half of these interactions. 

In eqs. (^-0) we explicitly wrote only a few representatives of the whole set of collision 
terms corresponding to the interactions listed in the table 1, the rest being denoted by 
ellipses. The c.c. terms that come from the evolution equation for njc are just the charge 
conjugates of the terms shown, and the term 

hdt-3H){nB-nBc) = ]:sB (9) 

O 

in eq. (|^) is due to the fast interactions of sphalerons, which change e^ number at 1/6 the 
rate at which they violate baryon number. (In eq. (P) we used the fact that s/s = —3H 
in the adiabatically expanding universe.) To see that this term is necessary, recall that in 
the absence of the AL = 2 lepton-violating interaction, we have the conserved quantity 
^B — Le = ^B — 2Le^ — Le^. If the AL = 2 scattering term in (|^) is zero, we see that this 
conservation law will be satisfied only if the ^sB term is present. Below we will express B 
in terms of Le^ and Le^. Lastly, let us remark that the second term in (|^) is particular to 
the model with electron lepton number violation and should be omitted when ^B — L^ is 
conserved. 

8 



Equations (§-0) constitute a set of nonlinear equations for the chemical potentials that 
appear in the distribution functions fi{p,fii). We can convert them into linear equations 
for the asymmetries by expanding the collision terms to first order in the /?/ij, which are 
proportional to the asymmetries through eq. (|^). Of course the resulting evolution equations 
are only valid for small chemical potentials, but this is a reasonable sacrifice because we are 
interested only in very small final asymmetries of the order of the present baryon asymmetry 
~ 10~^°. Even in the case of a large initial asymmetry the nonlinearity would persist only for 
a very short time, after which the system would be well-described by the linear equations. 
Then performing the expansions in chemical potentials and using equations (|l]) to rewrite 
the chiral fermion asymmetries in terms of /xq, equations (|§-|^) become 

-^ = - — (/^e« - f^e^ + /^O) x{Id + Ig + Y. h) 

^ = -2^-T,S(A„+^<.+2,„)W+jB (10) 

It should be noted that the chirality-changing part of the collision term in eq. (0) for Lg^ 
is proportional to /ig^ — /^e^ + yUo, clearly showing that Lg^ is driven towards the correct 
equilibrium value. These equations are not yet solvable because /iq and the sphaleron term 
|i? appearing on the r.h.s. of ([T0|) must be written in terms of the asymmetries Lg^ and Lg^. 
We will return to this issue shortly. 

The coefficients that multiply chemical potentials in ([10|) are defined as follows. The 
(inverse) decay process gives 

Id = 2/ndn, {27TyS\p^ - p2 - P3)\M{eL%R ^ h*)\'f^,f^,{l + f,,). (11) 

Scatterings involving the gauge bosons^ contribute a term 

Jg = 2 / n dn, i2n)^6\p, +P2-P3- Pa) x 
X {\M{e ,^en - h* A)\' f^J^,{l + f,,){l + /°J 

"^We thank Sonia Paban and Willy Fischler for originally pointing out this contribution to us. 



+ \MiAen - h*eL)\'f^^f%il - ^^ + fL) 

+ \Mihen - Ae,)\'f^,f,,il - 0(1 + /°j} , (12) 

where A refers to either of the vector bosons W or B. The scatterings involving fermions 
give 

If = 2 I n dn, i2TT)^6\p, +P2-P3- P4) X 

X {iMicL'^en) - /L7i?)r/°i/°2(l - /?3)(1 - /?4)+ 
+ ( IMihcR ^ fRCL)]' + IMifn'^en ^ Jl^cl)? ) x 

x[/°i/?2(l-/°3)(l-//4)]}- (13) 

Finally, the lepton number violating colhsion term is 



4 = 1 



j n dn,(27r)45^(pi+p2-P3-P4) 

[{Si\M{edL ^ h-h-)\^ + \M{eLU^ ^ h*h 



X ]^[ti(,\Jvi[eL(^L ^ IT' IT' )\^ + \JVl[eLUi,^ h'-h j|^)x 

x/ri/.°2(l + A°3)(l + /M) 

+ i\M{eLh'' ^ 4^/i-)|' + \M{eLh+ ^ i/A*)|' + |^(ei/i -^ z/A-)r)x 

X/eU°2(l-/.°3)(l + /.°4)}- (14) 

The distribution functions /°j appearing in ([TTl -|l4|) are the usual equilibrium functions de- 
fined in (H) for particle species a with zero chemical potential and momentum pi. 

There is a technical subtlety concerning the s-channel scattering contributions to the 
chirality-changing interactions. We have included the (inverse) decay contribution explicitly 
in our Boltzmann equations, which assumes that the Higgs particle is a thermalized member 
of the heat bath. However, the on-shell part of the Higgs exchange in the s-channel scattering 
physically corresponds to an inverse decay into a real Higgs boson followed immediately by 
a decay process. Thus the on-shell contribution in s-channel scattering introduces double 
counting of decay interactions and has to be subtracted away. Hence what we mean by 
the s-channel contributions in our Boltzmann equations is actually the on-shell subtracted 
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part. We will discuss the subtraction procedure in section 3, along with the computations 
of matrix elements and thermal averages. 

We still have the problem that in general /iq is not simply related to fie^ and /ie^. This 
would be the case for example if all the lepton numbers were violated by slow processes. Then 
four conditions would be needed to close the network of equations for the four softly- violated 
asymmetries 0]. But if we assume that only electron lepton number violation is slow, we 



can use the conservation laws for ^B — L^ and ^B — Lr (or the new chemical equilibrium 
conditions for /x^^ and /i^-^ if these are violated) to express the remaining unknown chemical 
potentials in ( p!0D in terms of /ig^ and /ig^. In particular we can parametrize /xq as 

fJ'O = ai^eR + bUeL (15) 

with a and b model-dependent constants, computable when the specific set of conservation 
laws and equilibrium conditions is given. 



Using (|]-|3) and ( |T5|) one then finds 



i? = ^(l-13a)L,,-l(l + 136)L,,. (16) 

By differentiating this with respect to time, we can write |-B in terms of Le^ and Lg^ in 
eqs. (p!0|), so that they close to give 



^en 



- r^z.((l + a)i^e^-(l-&)i^ej 

- 51 rAL,£(a£-^efl + fe£-^ei +Q), (17) 



dt 

dLe^ 5 + 13a dLg^ 

^r= ~ 13(1 +&) dt ,^^_^_^ 



where 



12 

/ 

Tal/ = 7f^I/lL,^- (18) 

In (p!7|), the coefficients a^ and 6^ are linear combinations of a and 6, and q are related to 
the initial values of conserved quantities. For example, if there is fast violation of muon and 
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tau lepton number, then /i^^ = ^r^ = ~/^0; as will be shown in section 5, and we obtain 

24 



24a 



b. 



<^ 13{l+b) ^ ]^3 

a^ = ar = \ae b^ = K = -be 

Ci = 0. (19) 



Nevertheless, in equations (ITTHTsD the parameters a, ai, b, and 6j are left unspecified in 
order to stress the fact of their model-dependence, notably their dependence on whatever 
boundary conditions are imposed by the fast processes. Since these appear multiplying 
the rates, the temperature at which the asymmetries reach their asymptotic values can be 
accurately determined only after these parameters are known. 

If there is no lepton number violation aside from sphaleron effects, that is if F^l = 0, the 
equation for Lg^ in (0) is actually redundant. In that case one would use the conservation 
of |-B — Lg along with eq. ([T6| ) to express L^^ in terms of Lg^ in the evolution equation for 
Lg^, which could then be solved for Lg^ directly; this is exactly what we will do in section 
4 below. If there does exist additional electron lepton number violation, \B — Lf. is not 



conserved and the more general form of (|T^ must be used. 



3 Rates for (inverse) decays and scatterings 

We are interested in the rate at which e^ number violation occurs in the early universe, to 
determine at what temperature the associated interactions come into thermal equilibrium, 
and in particular whether this can be close to the electroweak phase transition temperature. 
The dominant processes are usually the decays and inverse decays of Higgs bosons into e^ and 
the lepton doublet Lg shown in figure 1. As will be discussed in the next section, the answer 
to the question of whether an e^ asymmetry relaxes completely before sphalerons go out of 
equilibrium at the electroweak phase transition, is exponentially sensitive to the total rate 
of Ca violation. We have therefore undertaken to more carefully estimate this rate than was 
done in CDEO. We will find that for very small Higgs boson masses ?7i^, inverse decays by 
themsleves are slow enough so that e^ violation would not take place until after sphalerons 

12 



had largely disappeared, in which case the baryon asymmetry would survive despite the 
vanishing oi B — L. But before making this conclusion one should consider corrections to 
the rate of e^ violation due to two-body scattering processes. It will be shown that these 
actually dominate over the decays for small rriH, leading to a large enough rate of e^ violation 
that the results of CDEO are qualitatively unchanged. 

3.1 Decays 

For the (inverse) decay rate, we will make several refinements on the treatment given by 
CDEO: inclusion of the thermal masses of the leptons, use of the finite-temperature wave 
functions for the fermions, and retention of the Bose- Einstein statistical factor for the Higgs 
boson. The first two corrections are important in the limit of small m^ because they exac- 
erbate the the phase space suppression for the decay. 

The decay rate entering into the Boltzmann equation (|10|) is proportional to an integral 
Id, eq. ([Tl|) , in which appears the transition matrix element Ai. In terms of the electron 
Yukawa coupling h^, the squared matrix element for H -^ e^Cif^ is 

\M\^ = 2h%^-pe^ = 2hl{\pe^\\p,J -Pe^-Pe^), (20) 

where Pe^ and Pe^ are light-like four-momenta constructed from the corresponding three- 
momenta. These appear in place of the usual four-momenta because the fermion masses are 
from thermal effects which break Lorentz invariance but preserve chirality, so that one must 
still use massless spinors in constructing the amplitude [|ll|]. However to a good approxi- 
mation the dispersion relations for the particles appear to be massive. The thermal masses 
of the two chiralities of electrons and the Higgs, appropriate to the temperatures above the 



electroweak phase transition in which we are interested, are given by |11, 12 



ml^{T)/T^ = (3/ + ^'2)/32^0.044; 
ml^{T)/T^ = ^'78 = 0.017; 



■^We note that the coefficient of TOq in the expression for to^(T) is corrected here relative to that which 
was quoted in ref. Isl. 
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mi{T)/T' = 2D{l-T^/T')- 

D = {2mlr + ml + 2m'l + ml)/8v^. (21) 

Here rriQ is the zero-temperature Higgs boson mass, and g and g' are the respective SU{2)l 
and U(l)y couphngs evaluated at the running scale of 100 GeV which is at the lower limit 
of the temperatures we are interested in. Tc is the critical temperature for the electroweak 
transition, given at one loop by 

By substituting the lowest experimentally allowable values of the top quark mass rrit = 90 
GeV and rriQ = 60 GeV, we find that mu/T must be no smaller than 0.41 in the high temper- 
ature limit, where it becomes constant as a function of temperature. This also approximately 
coincides with the kinematic threshold for thermal Higgs boson decays into electrons. niu/T 
is the key parameter for the strength of the decays because Id depends entirely upon it, 
considering nte^ and me^^ to be known quantities. In Figure 2 we have plotted the thermal 
Higgs mass for various values of ttiq and mt, where it can be seen that niu varies between 
0.4 and 2 for mo between 60 and 1000 GeV, for relevant values of the top quark mass. 

The integral in eq. (pT]), which is proportional to the thermally averaged rate of decays, 
can be reduced to one dimension: 

_m„T^hl^'^ r°° e^mff/T f cosh.{ae^u + ^^u^ - 1) cosh(Q;e^M + ^^u^ ~ ^)\ 

TT^ Jl [e^rnn/T _ ly ycOsh(ae^M - 7a/m^ " 1) COsh(Q;e^M - '^^/u^ " 1)/ ' 

(23) 
where 

«eL = {ml + ml^ - mlj/AniHT; 
"efl = {ml + ml^ - mlJ/ArriHT; 
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\'/\ml,ml^,mlj/Am„T; 



\{x,y,z) = {x-y-zf-4yz. (24) 

In CDEO, Id was approximated by {herriHT In 2y/{16TT^). In Table 2 we compare the nu- 
merically integrated value with this approximation for some values of m^/T. 
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m^/T 


Id,cdeo/K 


Id/K 


hop/K 


Ig/K 


0.4 


1.6x10-^ 


2.5x10-5 


4.2 X 10-5 


3.2x10-^ 


0.6 


3.5x10-^ 


2.0x10-^ 


4.2 X 10-5 




0.8 


6.2x10-^ 


4.1x10-^ 


4.1 X 10-5 




1.0 


1.0x10"^ 


6.4x10-^ 


4.0 X 10-5 




1.2 


1.4xl0~3 


8.5x10-^ 


3.7 X 10-5 




1.4 


1.9x10"^ 


1.1x10-3 


3.4 X 10-5 




1.6 


2.5x10-3 


1.2x10-3 


3.2 X 10-5 




1.8 


3.1x10-3 


1.4x10-3 


2.9 X 10-5 




2.0 


3.9x10-3 


1.5x10-3 


2.7 X 10-5 





hlT' 



Table 2: The quantities /„ that appear in the Boltzmann equation, divided by K 
The subscript CDEO refers to the approximations for the decay rate made in reference |^, 
while D, top, and G refer respectively to the exact decay rate and the scattering rates for 
top quarks and gauge bosons to produce e^. 



3.2 Scatterings producing Cr 

One might expect that two-body scatterings which change e^ number are slower than decays 
because of phase space and coupling constant suppression. However in the region of small 
thermal Higgs boson masses the phase space for the decay is strongly suppressed. In this 
regime the production of electrons by top quark scattering (Figure 3) or Higgs/gauge boson 
scattering (Figure 4) may actually dominate over the decays. 

In the Boltzmann equation, the thermally averaged scattering rates for processes which 
change e^ number appear through the integrals // and Iq defined in eqs. (|T^|-^). We may 
ignore all the fermion scattering processes except for / = top since the top quark has the 
dominant Yukawa coupling. The channels which produce an e^ in the final state are listed 
in Table 1, and are categorized here as top-quark or gauge/Higgs boson scatterings. We 
use the Maxwell- Boltzmann approximation for the initial state particles, and neglect the 
distribution functions compared to unity for the final state particles. The errors introduced 
by this procedure are irrelevant in the regime of large m^/T where the scatterings are small 
corrections to the decay rate; in the small mu/T regime, although the scatterings dominate. 
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they are still sufficiently fast for our purposes that the anticipated 20% corrections to the 
Maxwell-Boltzmann approximation would not change any of our subsequent conclusions. 
Then the quantities li can be expressed as an integral over the Mandelstam variable s 
involving the corresponding cross sections and the modified Bessel function Ki. 

Ir = ^^E/,,, -^^\{s,mlml)K,{^s/T)a^{s) (25) 

The sum is over all possible channels c, each of which has a threshold determined by the 
thermal masses of the initial or final state particles, M^ = max(mi + 7712,^3 + 1714). The 
sum includes a factor of 2 for isospin (since all the processes we consider involve doublets) 
and if there are gauge bosons, a factor of 3 for the polarizations. Note that there is no 
corresponding factor for the fermion spins because they all have a definite helicity. The sum 
also runs over the different crossings that leave e^ in the final state. The choice of whether 
Cfl is in the initial or final state is arbitrary-either way gives the same final answer. 

In performing the thermal averages of the cross sections, we omitted corrections of order 
m^/s, where m is a thermal mass, since typically m < T and s ~ 20T^. Thus we can set the 
threshold to zero in the integration of eq. (^) and approximate X{s, ml, m^) by s^. (The 
same neglect of thermal masses would not have been appropriate for accurately determining 
the decay rate since the phase space there is determined by the scale 771^ rather than s.) In 
this approximation, the analytic expressions for the top quark scattering cross sections are 

and for gauge/Higgs boson scattering 

E/ \ e /2i 2i 2 21 /r\'-T\ 

, ^^(^) = ;^ l^-e. In ^ + m,^ In -^ - -m,^ - -m,^ j , (27) 

where the thermal masses of eq. (|2T|) have been used and F is the width of a Higgs boson 
with mass y/s: r(s) = ^^V^- Some care must be taken to obtain the subleading terms 
correctly. The numerical values of the above contributions to the Boltzmann equation for 
the Cjj asymmetry are also shown in Table 2. We have assumed a top quark Yukawa coupling 
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corresponding to a mass of 90 GeV in computing Itop- Although this may be unreahstically 
small, we wanted to examine the regime of small itlh/T which is most favorable to the 
possibility of keeping e^-changing interactions out of equilibrium below the freeze-out of 
sphalerons, and eq. (^) shows that this occurs for the smallest values oimt and itiq. However, 
Table 2 shows that the contributions from scattering dominate those of decays in the low-m^ 
region, so that one cannot in fact evade the washout of an e^ asymmetry even by going to 
small masses. More will be said about this in section 4. 

3.3 Subtraction of Resonances 

As has been mentioned above, the s-channel scattering tflti'^ -^ GrGl'^ contains an on-shell 
part which must be subtracted to avoid double-counting the contribution from the inverse 



decays. This necessity has been mentioned in references [^ and ||1J], although the details 
of how to carry it out were not explained. We have considered two methods that seem 
reasonable and are in good quantitative agreement with each other. Writing the finite- 
width, squared s-channel propagator for the Higgs as P(s; V) = ((s — m|)'^ + sF^)"^, we can 
summarize the different methods by making the following modifications to the propagator, 
which we refer to respectively as the principal value (PV) and delta function subtraction 
(SDF) methods: 



PV: 


P{s 


r) 


^P{s 


+ ie; 


0); 






SDF: 


Pis 


r) 


^P{s 


n- 


■>"- 


2 \ 


(28) 


hod gives /^„ 


X-2 


dx = 


-a-i 


— 6^^, even when the double 


pole is 



straddled by —a and b. The SDF method attempts to isolate the off-shell part by subtracting 
what the scattering rate "would have been" if the resonance were infinitesimally narrow. 
Although these alternatives at first look quite different from each other, by writing 7r(5(s— ?7i^) 
as e/((s — wi^)^ + e^) one can show that they are actually equal to each other in the limit as 
F — >• 0; writing x = {s — m|), and e' = eF, the difference is 

versus T-ir—?TTr-Z?TT^.- (29) 



(x2 + e'2)2 (a;2 + e'2)2 + x2F2 
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m^/T 


PV 


SDF 


0.4 


1.52 X 10-5 


1.52 X 10-5 


0.8 


1.43 X 10-5 


1.42 X 10-5 


1.2 


1.04 X 10-5 


1.04 X 10-5 


1.6 


5.24 X 10-^ 


5.33 X 10-6 


2.0 


3.95 X 10-^ 


3.58 X 10-7 


2.4 


-3.53 X 10-6 


-3.53 X 10-6 


2.8 


-6.15 X 10-6 


-6.16 X 10-6 


3.2 


-7.14 X 10-6 


-7.55 X 10-6 


3.6 


-7.70 X 10-6 


-7.67 X 10-6 



Table 3: Contributions to hop/K (see table 2) due to s-channel Higgs boson scattering. We 
compare the principal value prescription used in table 2 to the delta function subtraction 
method. 



In fact this is as good a correspondence as one should hope for since the whole concept of 
asymptotic states which have thermal distributions breaks down if they become very broad 
resonances. Thus two subtraction schemes should be considered equally good if they differ 
only by terms of order F^. We compare the PV and SDF methods in Table 3, of which the 
former was used in the top-quark contributions shown in Table 2. 

The negative values for mfj/T > 2.0 may be alarming at first, since they are supposed 
to be the off-shell part of a scattering rate. Although our results are completely insensitive 
to the negative contributions because of their relative smallness compared to the decay rate, 
it is interesting to speculate on their meaning. We take them as being suggestive that 
the decay rate Id gives an overestimate of the Higgs-mediated processes at large values of 
ttih/T. In this regime, due to strong couplings to itself or the top quark, the Higgs boson 
is starting to become a rather broad resonance where it may not be sensible to consider the 
decays; rather one should include the scatterings only, without subtracting the resonance 
from the s-channel. These give a smaller contribution than the decays because the area 
under the resonance curve gets diminished by the Boltzmann suppression in the function 
i^i(yi/T) ~ exp{—^/s/T) (eq. (p5|)). But since the off-shell scatterings are completely 



negligible compared to the decays in the large m^j/T regime, the question of whether to take 
the negative cross-section seriously is moot-its contribution is negligible. 

3.4 AL = 2 Scatterings 

In a later section we will also be interested in how AL = 2 lepton number violating inter- 
actions come into equilibrium, in connection with baryon number erasure. The operator we 
consider, which gives rise to neutrino masses, has dimension five: 

^-J(L,if)(l + 75)(//^^-)+h.c., (30) 

ij 

where rriij is the neutrino mass matrix in the fiavor basis, v is the 246 GeV VEV of the Higgs 
field H, and Li are the lepton doublets. 

The sum over channels of the cross sections for the processes e^^z/^ —^ h*h^ and e^ii —>■ 
h~h~ must be properly weighted by the factors Si shown in the Boltzmann equation (0) for 
Cl, because of the complication of having identical particles in some of the initial or final 
states. We find that 

. "-^(^) = 4^^ 1 6K,P, £ = /.,r. ^31) 

This is what we must substitute into eq. (|25| ) to obtain the quantities Ial/, which appear 
in eq. ( pi]) for the rate of lepton violation. The various channels respectively contribute the 
relative weights (2, 2,4, 1, l)|meep for 6^6^ —>■ h~h~, e^Ve —>■ h~h*, e^h^ —>■ e^fh' , e^h^ -^ 
Ve'^h*, e^h -^ Ve^h~ and (4, 2,4, 1, l)|me£p for e^l^ -^ h~h~, e^z/^ -^ h^h*, e^h^ —>■ i^^h', 
ej^h~^ -^ i>£^h*,ei^h -^ i'e'^h~ with i = fi,T. In the approximation of ignoring masses compared 



to s, the integral in eq. (|25D can be done exactly, 



T^6 

^AL,£ = ^- y] (^c,e- (32) 



" c 



This simplification occurs because of the constant cross section. Of course at very high 
temperatures of order the mass M of the heavy neutrino which must have induced the 
operator (^), the cross section will exhibit a form factor behaving as M"^ / {T"^ + M^) , coming 
from the propagator. Thus we are assuming that M is greater than the temperatures of 
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interest when we use eq. (p2D. It will be shown later that as long as M exceeds 10 TeV, our 
neglect of the propagator is justified. 

4 The case of conserved B — L 

The importance of the weak electron chirality-changing interactons in the context of pre- 
serving the baryon asymmetry was first discussed by the authors of ref. [Q (CDEO). They 
considered whether in a universe with primordial B — L = 0, the right-handed electrons 
could stay out of equilibrium long enough so that the primordial baryon asymmetry, kept 
nonzero by the asymmetry in the e/j species, does not get washed out completely before 
the sphaleron interactions fall out of equilibrium. While their result was negative for the 
preservation of the baryon asymmetry, they noted that the outcome of the computation is 
exponentially sensitive to the strength of the chirality-changing interactions. Therefore even 
a rather small change in the interaction rates could alter the conclusion. 

Indeed, the derivation of CDEO contained a number of inaccuracies with the potential to 
change the outcome significantly. In particular, their neglect of the final state particle masses 
leads to an overestimate of the rate of Higgs decays, and therefore to an underestimate of 
the final asymmetry. As we have seen in the section 3 above, the finite mass effects are very 
important for small (thermal) Higgs boson masses. On the other hand, CDEO neglected 
all scattering processes, which we found to be of particular importantance when the decay 
process is suppressed by the phase space effects. Armed with our accurate expressions for 
all tree level chirality changing rates computed in the previous section, as well as with the 
carefully formulated evolution equations derived in section 2, we can now re-examine and 
answer definitively the question CDEO posed. 

As explained in section 3, the evolution equations for the asymmetries depend on the 
particular set of equilibrium conditions and conservation laws characterizing the system. 
Here the set of boundary conditions is given by the equilibrium equations (|l]-0) supplemented 
with the additional constraint of vanishing B — L. In fact we only assume that primordially 
{B — L)p = 0; the vanishing of -B — L at any lower temperature is guaranteed by the first 
three of the conservation laws (^, which also imply the conservation oi B — L. In terms of 
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the chemical potentials we then find the constraint 

B - L = 12/i„^ - 3 5Z A^^L - 2/iei - f^cR + 2Aio = 0. (33) 

£ 

From eqs. (0-0) and (|33|) one can readily find that 

5 



yWo = 


- ^53(/ie« /ie. 




44 ^ 


Ub = 


- ^53 (/^e, ^e. 



(34) 

Note that if we further assumed the cr population is in equilibrium, i.e. fie^ — /^e^ = — /^o, 
then only the trivial solution would exist: fiQ = B = L = 0, in agreement with the standard 
analysis of ref . [HI . 



Our goal is to follow the evolution of the e^ asymmetry using eq. (|iy), where from (|3^ ) 
we see that the constants a, b are given by a = —b = 5/153. However, we need not consider 
coupled equations for Lg^ and Le^ because the latter is determined through the conservation 
law |i? — Lg = Ce^p- Eliminating B from this equation using (|3^ for fis oc B, we find that 

i,^ = _115i,, _ 15£c„. (35) 

962 ^ 962 '^ ^ ' 

Inserting this into eq. ([T^) then gives the simple first order equation for Le^ 

^ = -nLen + \Ce,p). (36) 

where 

711 17.7 

r = -^Trl = ^Id + Ig + hop). (37) 

The quantities Ix/T^ are the thermally averaged rates computed in section 2. Eq. ( PBj ) is 
easily integrated to yield a general solution for Lg^ as a function of time that is completely 
specified in terms of primordial quantities: 

LeAt) = -^a,p + {Le,,p + ^a,p)e"^'o'^*'^. (38) 

Interestingly, the right handed electron asymmetry does not tend to zero, but instead 
relaxes to a fixed value — |Ce,p which may be even larger in absolute value than the primordial 
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value of Le^. This is contrary to the assertion of CDEO that Le^ -^ 0. What happens 
is that at the fixed point the right and left chiral asymmetries become equal under the 
additional constraint of (|-B — Le)-conservation. This results in erasing the baryon and 



lepton asymmetries, B, L ^ 0, as can be inferred from eq. q^ It is useful to write down an 
expression for B directly as a function of time 

i?(t) = 5e,e-^'>*'^, (39) 

where the equilibrium value of B much above T* is given by 

66 

B,, = —{a,, + 3L,,,,). (40) 

As expected, the value of B decreases exponentially in time, simultaneously with the equi- 
libration of the e/j population. It is conceivable that the initial value for B^g could be as 
large as B'^'^ ~ 10~^ in some grand unified theories, although typically it will be suppressed 
by small coupling constants and CP-violating phases. We will allow B'^'^ to vary in what 
follows. 

Equation (|39|) is our main result in this section. Although we will compute the rate F and 



its time integral for (|39D numerically, to get a feeling for the relevant scales it is useful to find 
an analytic approximation for -Bgnai- We can do this by employing linear fits to the values 
of the functions Ix given in table 2, and integrating the ensuing fit functions analytically. 
Using a simple linear regression formula, we find 

Trl - 3.8 X 10-'Thlf{m^{T)/T), (41) 

where the fit function 

f{x) = (-1.1 + 3.0x) + 1.0 + /ij2(0.6 - O.lx) (42) 

has been normalized so that the middle term, due to gauge/Higgs scattering, is unity. Be- 
cause the number of degrees of freedom is constant, g^, = 106.75, entropy is conserved and we 
can use the approximation T/T = —H. Then, neglecting the threshold effect in the decays. 



we obtain the result 

'■*- ^100GeV\ 



T' dtr ~ 350 ( TfA-) [(-1-1 + 2.40;^) + 1-0 + /i?(0.6 - O.OQx^)] , (43) 
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where we have defined the asymptotic high-temperature hmit of the Higgs boson thermal 
mass as 

hm m^(T)/r = V2D = xh (44) 

T — ^oo 

The approximation (^31) treats the decay part too roughly for small thermal Higgs masses, 
and therefore cannot be used to obtain figure 4a. However it is a very good estimate for the 
total integrated rate over the whole interesting range of x// <,2. It should be noted that for a 
large range of thermal Higgs masses, xh^O.9, the scattering contribution to the integrated 
rate dominates over the decay part. 

In figure 5 we show the result of our computation for the final baryon asymmetry in units 
of the initial asymmetry (^), as a function of the physical Higgs and top quark masses. For 
comparision to the work of CDEO, we show in figure 5a what the outcome would have been, 
had we included only the corrected decay contribution to the total rate. The curves are 
contours of constant -Bfinai- The initial upward bending of the curves as a function of the 
zero-temperature Higgs boson mass rriQ occurs because the endpoint of the integration, the 
electroweak phase transition temperature Tc, increases with mo. The resulting decrease in 
the integrated rate can be compensated by increasing nit. As tuq continues to increase, the 
thermal Higgs mass starts to show its dependence on mo (cf. eq. ^l]) leading to a significant 
increase in the rate as a function of mo. Simultaneously the phase transition temperature 
reaches a plateau and therefore the contours turn down at large values of m,o- 

It is noteworthy, and in contrast to the result of CDEO, that their mechanism could 
have appeared to be successful in protecting the baryon asymmetry for a small region of 
experimentally allowed parameter space, had they considered only the more exact expression 
for the Higgs boson decay rate. The main reason for this would-be positive result is the 
suppression of the decay rate for small values of xh, where the phase space volume is small. 
When the scattering contributions are included however, the conclusion changes drastically. 
We show the result of the numerical integration in figure 5b. In particular, because of 
the gauge scatterings the final baryon asymmetry is completely negligible over the entire 
parameter space. Thus the scatterings play a decisive role in the most relevant region of 
parameters for the problem. In the regime where the decay contribution dominates over the 
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scattering rate, -Bfinai is even more suppressed than in the region of parameters shown in 
figure 5b, which extends only up to xh ~ 0.55. 

Our conclusion is qualitatively the same as that of CDEO: it is not possible to protect 
the baryon asymmetry of a 5 — L = universe by temporarily storing it in the asymmetry of 
the cr species. However, the way in which we reached this conclusion differed from theirs in 
important respects and moreover, the robustness of our result erases any doubts expressed 
by CDEO, that some subtle thermal effects might alter the outcome. 

5 Examples when B — L is not conserved 

Although we have seen that an e^ asymmetry cannot save the baryon asymmetry from 
sphaleron effects when B — L = 0, it can have a dramatic effect when B — L ^ 0, as 
one might have in SO(IO) grand unified theories where B — L gets broken at the GUT 
scale. In this section we will explore the consequences of an e^ asymmetry in conjunction 
with sphalerons and other B- and L-violating interactions that hitherto were believed to 
completely erase B. 

It is expected that any global symmetries in nature are only approximate; they might be 
broken by exotic gravitational means such as microscopic black holes or wormholes, or more 
mundane reasons like the fact that one can construct renormalizable B- and L-violating 
interactions in supersymmetry using the squark and slepton fields. Such effects may be 
described by nonrenormalizable operators in the low energy effective theory. Assuming that 
a primordial baryon asymmetry survives allows us to constrain the strength of these new 
interactions. These constraints were first applied to the AL = 2 operator (^) by Fukugita 
and Yanagida 0]. They were generalized to other operators in refs. [^, ^. 

Quite generally, an operator of dimension 4 + n would have a dimensionful coefficient 
of the form M~", and would mediate a process at high temperatures whose rate scales like 
r ~ 2^2n+ij^^-2n_ This rate is to be compared to the Hubble expansion rate ~ T'^/Mpi. 
Insisting that the interaction be out of equilibrium below the temperature Tq at which 
sphalerons came into equilibrium (10^^ GeV), or when the baryon asymmetry was formed, 
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whichever was smaller, leads to a lower bound on the scale of new physics: 

M>r„(^)"". ,45) 



In supersymmetric models it was argued ||T5| that due to additional anomalies which can 
temporarily protect the asymmetry (until the effects of supersymmetry breaking kick in), the 
maximum temperature should be at To ~ 10^ GeV rather than Tq ~ 10^^ GeV. Interestingly, 
in the context of inflation the maximum temperature should not surpass the thermalization 
temperature [0, which could be as low as ~ 10^ GeV IT^. We will argue that actually 
none of these values of Tq are the correct temperature to use in setting this bound whenever 
an Cfl asymmetry is present, as it will prevent the baryon asymmetry from relaxing to zero, 
essentially because of the charge neutrality of the universe. Rather, one should use the 
temperature T* at which e^-violation comes into equilibrium and allows the e^ asymmetry 
to disappear. In this section we will compute T, ^ 1 TeV. Hence constraints like (|45|) might 
be weakened by as many as nine orders of magnitude, by standard model physics which has 
been known for well over a decade! Note, however, that our arguments here do not affect 
the limits imposed on renormalizable B- or L-violating operators. For these, the strongest 
bound always comes from applying the limit at Tc (when sphalerons disappear) [^. 

Before computing T^,, we will follow the evolution of the baryon and e^ asymmetries 
more exactly, for a particular example of a nonrenormalizable operator, that which would 
be responsible for the electron neutrino mass. An example employing explicit i?-violation 
will also be given. 

5.1 Neutrino masses 

The best known example oi B — L violation is the possibility of neutrino Majorana masses, 
which would be generated by the dimension five operator, ^mijV~'^{LiH){LjH). To avoid 
an erasure of the baryon asymmetry, upper limits on mjj, the neutrino mass matrix elements, 
must be satisfied P, ||, |18[ . The flavor-dependence of the operator will be important since 
any approximately conservered quantities (such as ^B — L{) would protect the asymmetry 
|T8| , |T^ . If for all flavors the operator was so weak as never to have been in equilibrium, and 
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B — L = 0, then as shown in the previous section, the baryon asymmetry would be erased, up 



to small mass effects I^O], |T9| . If, on the other hand, one or two of the lepton generations were 
violated by this operator and the remaining generations were out of equilibrium, then even 
ii B = L = (or more generally, B — L = 0) a baryon number would be generated provided 
some initial lepton flavor asymmetry existed p|. But if all three flavors were violated down to 
low enough temperatures, the baryon asymmetry would be erased. How low is low enough? 
The answer, which we will denote by T*, depends crucially on the decouphng temperature 
of the chirality-fiipping electron interactions. In this section we will determine T*, as well as 
the size of the neutrino mass matrix elements which would be compatible with preserving a 
primordial baryon asymmetry. 

Supposing that the masses of u^^ and u^ were sufficiently large (several keV) that the 
separate lepton flavors L^ and Lj. were strongly violated during the sphaleron epoch, we 
would insist that L^ is not similarly violated, if we want to keep sphalerons and neutrino 
mass effects from erasing the baryon asymmetry. Since we are working above the temperature 
at which the Higgs gets its VEV, however, the effect of interest will be scatterings of leptons 
into bosons. 

By assumption, the AL = 2 scattering processes LiLj -^ HH will be in equilibrium if i 
and j are either /i or r, which means that in place of the conservation laws for |S — L^ and 
^B — Lr in eq. (^, we will have two nonstandard equilibrium conditions 

/^Mt = /^Ti = -/^o (46) 

involving the Higgs chemical potential. As in the previous section, the new constraints allow 



us to solve for the parameters a and b in eq. (|T^). We find that 
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165^^« ~ n^ 



/is = T^/^efl-T7/iei- (47) 



In addition to the equilibrium conditions (|46|) for the mu and tau species, the new interac- 
tion provides a source of electron-type lepton number violation which goes out of equilibrium 
at a temperature depending on the size of the neutrino mass matrix elements rriei- The effect 
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is parametrized by the integrals /al,j in the Bohzmann equation (|T0|), calculated in section 
3. The evolution equations for the two electron chiralities then become 
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(49) 

Using the results of section 3, we have integrated these equations starting at high temper- 
atures T ^ 100 TeV where the e^- violating interactions have not yet reached equilibrium. 
We consider a range of initial conditions for Le^ , including the largest reasonably imaginable 
asymmetry in e^, which would occur if there were only e^" and no e^"*"; then dividing ng^ 
by the total entropy gives Le^(O) = 10^^. The initial value for Le^ is given by solving the 
chemical potential equations subject to the constraint that the Lg-violating interactions are 
in equilibrium: /ig^ = — /io, as in eq. (|4B|). We find that Lg^ = — ^Lg^ and B = \Lf^^ initially 
(see section 5.1, below). The resulting value of the baryon asymmetry is then determined 
by the final lepton asymmetries using eq. (^T]) : B = ^-Z^ej, — if-^e^,- 

In figures 6a and 6b we show the evolution of the asymmetries Lg^, Lg^ and i? as a 
function of the temperature for rrit = itlh = 100 GeV. In figure 6a we have assumed that 
the electron lepton number violation has gone out of equilibrium before the e^-e^ transitions 
come into equilibrium. Initially all asymmetries have their high temperature equilibrium 
values. As the right handed electrons come into equilibrium, the asymmetries relax to their 
low-temperature equilibrium values, which are of the same order as the asymmetries above 
the transition region and are calculable in terms of the primordial Lg^ asymmetry. In figure 
6b we assume that the neutrino mass parameter 

(mgj)^ = -|mge|^ + |mg^|^ + \mer\^ (50) 



"The relative coefficient of |?neeP, which was wrong in ref. pj, has been corrected here. 
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takes the value (rriej) = 5 keV. Now the asymmetries fall exponentially within the region 
where the lepton number violation and the Cn-e^ transitions (with rate Fe^^e^) are both in 
equilibrium. A new equilibrium is reached above the electroweak phase transition tempera- 
ture, Tc ~ 154 GeV, as the lepton number violating interactions fall out of equilibrium. The 
lesson to be learned from figure 6b is that, in order to erase the asymmetries, the tempera- 
ture must fall somewhat below that at which Fe^^eiZ-f^ = 1, as is often assumed in simple 
estimates. 

To obtain a conservative bound on {ruei), we must assume that the lepton number vio- 
lation is strong enough to decrease the largest imaginable initial baryon asymmetry down 
to the smallest value of the present baryon asymmetry -Bmm — 4 x 10^^^ inferred from big 



bang nucleosynthesis [21]. For the former value we take -Bjnitiai = 2 x 10 , corresponding 



to Lg^^p = 10~^. Moreover, we will assume the least favorable values rrit = 90 GeV and 
tuh = 60 GeV for the unknown masses of the top quark and the Higgs boson. In this way 
one obtains the bound 

Ki) <20 keV, (51) 

an order of magnitude less stringent than the rough estimate given in our letter 0. This 
result can also be inferred from figure 7, where we show the results of a more complete 
computation of the bound as a function of the initial asymmetry. We also display the 
dispersion in the bound that results from the above-mentioned uncertainty in the present 
baryon asymmetry and the top quark and Higgs boson masses. The shaded area corresponds 
to the intervals 90 GeV < nit < 300 GeV, 60 GeV < m^ < 1 TeV and 4 x 10"^^ < 5finai < 
6 X 10-". 

A simple analytic estimate for the bound (pi]) can be obtained by requiring that the rate 
of lepton number violation, or more generally, any other baryon or lepton number violation 
under consideration, becomes comparable to the expansion rate after the time t^, defined by 

/** dt Tlr = ln(Sinitiai/4 X 10-") ~ 18, (52) 
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using our assumption that -Binitiai = 2 x 10 '^. To solve equation ( |5^ ) for the correspond- 
ing temperature T^, we use the linear fit for Tlr derived in section 4. After analytically 
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integrating the fit function and expanding in tlie small parameter T^jT^ we find that 

T, ^ 1.3/(x^) TeV, (53) 



where f{xu) is as in eq. (|43|) . Again assuming that mt = 90 GeV and tuh = 60 GeV, we 
get xh — 0.41 and then T,. ^ 1 TeV. Note that while we have so far loosely spoken of T* 
as the temperature at which e^-e^ transitions come into equilibrium, it is actually given by 
the considerably smaller temperature (p3D, which is relevant for obtaining the bounds on 
various baryon or lepton number violating operators. Let us finally note that using (|5^) in 
the rough condition for freezeout of the lepton-violating interactions, 

Tal 



H 



T=n 



1, (54) 



gives a bound (mej)^21 keV, in good agreement with the accurate numerical result (pT]). 

5.2 Baryon number violating operators 

Apart from neutrino masses, there is a host of possible {B — L)-violating interactions that 
have explicit baryon number violation. To further illustrate the relevance of right-handed 
electrons for preserving a baryon asymmetry, we will show what happens when several such 
operators are in equilibrium during some epoch. 

It is convenient to rewrite the sphaleron and electric charge constraints on the chemical 
potentials in the form 



f-^- ^ 


-- 2^e + -{L, + Lr) + -L,^; 


l"^> ^ 


-- -lLe-^{L, + U) + h,^, 



(55) 

where H is the Higgs boson asymmetry, L^ = 2Le^ + Lg^ is the asymmetry for total electron- 
type lepton number, and similarly for the other flavors. We have noted that in the minimal 
standard model, there are three conserved quantities ^B — Li and one which is approximately 
conserved at high temperatures, Lg^. In the absence of new interactions we would solve for 
the baryon and Higgs chemical potentials in terms of the conserved quantities. 
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When sufficiently fast interactions are introduced which violate any of the conservation 
laws, eqs. (|55| ) must be supplemented by the new chemical equilibrium conditions associated 
with these interactions. Of course we do not want to break all four of the symmetries, since 
then the solution to the system of equations would be the vanishing of all asymmetries. We 
do wish to break the particular linear combination B — L, however, since only in this case do 
we get anything new from our observations concerning right-handed electrons. If 5 — L is 
conserved during the epoch when sphalerons are in equilibrium, yet has a nonzero value, this 
by itself is enough to insure a nonvanishing baryon asymmetry, regardless of the electron |^ . 

For example, consider the AB = 2 operator 

M-\u^d%){d^dl){u^dl). (56) 

It conserves the linear combination 

2Le -L^-Lr = C. (57) 

(Although it conserves both Lg — L^ and L^ — Lj. separately, only the sum of these is relevant 
since L^ and Lr enter eqs. (|55| ) in the combination L^ + Lr.) It also leads to the equihbrium 
constraint 

2/i„^ + 4/^rf^ = ^ B = 4H, (58) 

where we used eq. (|I]) to eliminate the chemical potentials of the right-handed quarks in 
favor of fiqj^ and the relation ^b = 12/ig^. It is now possible to solve for B and the total 
lepton number L using eqs. (|55| , |57| - p8D , and in particular their difference^] 



3 

68' 



5-L = — (C-9L,^,p), (59) 



where we subscripted Le^ to emphasize that it retains its primordial value until T*, the 
temperature at which the rate of e^-violation starts to exceed the expansion rate of the 
universe. Let us suppose that the new interaction (|56D goes out of equilibrium at some 
temperature T' below the baryogenesis scale but above T*. Then for T' > T > T^, the 
^This is similar to the result in ref. but written in terms of the conserved quantities C and Lg^. 
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asymmetry B — L maintains its equilibrium value ([59|), since nothing occurs that would 
change it. As the temperature falls below T*, B — L continues to be conserved, even though 
Le^ itself changes. At this point, the complete equilibrium conditions assumed by Harvey 
and Turner prevail, and we can use their results to determine the final B in terms oi B — L. 

This result is remarkable in two respects. First, B — L in (|59|) is completely independent 
of its initial value, at temperatures before the AB = 2 operator came into equilibrium. One 
possibility is that it was initially zero, in which case the new operator actually worked in 
conjunction with sphalerons to create the baryon asymmetry, which would have been zero 
in its absence. So we have here a new way to convert a primordial lepton asymmetry into 
baryons. Secondly, even though this mechanism depends on the conservation of e^ number, 
it does not require that Lg^ be nonzero. It was enough to have an asymmetry between the 
different flavors of leptons in this example. 

In general we can violate all three of the quantities |i? — Li and still obtain a nontrivial 
solution for B. For our second example, we take the operators 

M-5(M«u/)(u«£,,/)(L,r2Lfc'=), (60) 

where the indices i,j,k are generational and T2 is the Pauli matrix. Although there are 
numerous possible combinations of i,j,k, consider at first any three which break all three 
symmetries. These give relations among the chemical potentials which as before can be 
expressed entirely in terms of left-handed fields using the constraints (|ID from Higgs boson 
interactions to eliminate the right-handed fields. After so doing, the new relations are 

3/Xg^ = -/i^L = -f^TL = /^o, (61) 

which translate directly to 

B = -L^ = -L^ = AH (62) 

in terms of the asymmetries. We can solve these together with eqs. ( pSj ) for all of the 
asymmetries in terms of Le^ : 



1 
20' 



H = T^^ea, 
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B — —J^fj, — —i-^T — —i-'en'i 

9 

Le = Y^^en- (63) 

The final baryon asymmetry will again be determined hj B — L, whose equilibrium value 
above T* is 

B-L = -^L,^,,. (64) 

as first shown in [^. 

If the primordial asymmetry in e^ was sufficiently large, the new B- and L-violating 
interactions (|56| , |60D were harmless for the baryon asymmetry, so long as they fell out of 
equilibrium above T^. As outlined in the introduction to this section, this leads to upper 
bounds on the masses in their coefficients [^. Let us examine the bound for the dimension 
nine operators M~^qqqqqq and M~^qqqlll we have just discussed. For a particular scattering 
channel, say for a 3 ^ 3 process, we saw in section 2 that the relevant rate is typically given 
by r ~ 12//T3, where 

/ = /ndn.(27r)^5(pi +P2 +P3 -P4 -P5 -P6)|-Mre~'3(^i+^^+^^) (65) 

and we can estimate the matrix element as 

i 

There are 6-5-4 different channels for the 3 — > 3 processes and 6-4 for the 2-^4 processes. 
Putting everything together and evaluating the integrals, we ffnd a rate of 

r = 288(2/7r)^^T(^-j , (67) 

which, if we insist it be less than the Hubble rate at temperature Tq, gives the lower limit 
on M of 



M > O.TTo ( ^ ) , (68) 
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in remarkably good agreement with the dimensional analysis estimate (^S]). Previous limits 
on this operator [Q ranged from 10^ — 10^^ GeV for Tq between 10^ — 10^^ GeV. We have 
shown that the appropriate temperature is Tq = T^, ~ 1 TeV. The difference for the bound 
on the scale of new physics is marked: 

and hence the limit is much closer to the weaker of the bounds given in |Q . Limits on other 
nonrenormalizable operators are similarly weakened. 

The reader may wonder how general our results are. Did we need exactly three of the 
operators (^) in equilibrium, and did it matter which three? Intuitively we expect that the 
answer to both questions is no as long as all three symmetries |-B — Li are broken, since 
the only thing keeping the chemical potentials from vanishing is the necessity of maintaining 
charge neutrality of the universe by canceling the charge of the conserved e^ asymmetry. 
Therefore it should not matter how many or which kind of operators do the job of breaking 
3-B — Li. 

It is easy to prove that having other operators in addition to all those we have already 
considered has no effect on the equilibrium conditions (0). Imagine some operator which is 
perfectly arbitrary except that it conserves e^ number. It can be characterized by a set of 
integers np, the net number of fields of type F, defined to be the number of fields F minus 
the number of charge-conjugated fields F'^. If it is in equilibrium, it gives a condition on the 
chemical potentials, 

Y,np^p = Q. (70) 

F 

But the coefficients np must be such that the total hypercharge of the operator vanishes, 
because of gauge invariance: 

Y.yFnp = Q. (71) 

F 

Now note that if the solution for the chemical potentials takes the form ^p = cYp |p2|| , 
then these equations are consistent with each other. In this case, eq. ( [70D is automatically 
satisfied, so it is not really a new condition, as we set out to demonstrate. It is easy to 
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verify from eqs. ( |63| ) that the chemical potential of each particle is indeed proportional to 
its hypercharge. Therefore we can add as many new operators as we want: as long as they 
conserve e^ number, they lead to the same equilibrium conditions, which can be written as 



1 
20 



l^A = T^^A/^efl- (72) 



6 Summary and Conclusions 

Although the smallness of the electron Yukawa coupling is not completely understood, its 
smallness is an experimental fact. In the absence of new physics beyond the standard model 
and still below the GUT scale for which the right-handed electron is not a singlet, e^ remains 
an essentially decoupled state for a good part of the early evolution of the Universe. Indeed 
we have accurately determined the decoupling temperature of this state to be T^ = 1.3f{xH) 
TeV. 

Our interest in the e^ decoupling temperature is directly related to the preservation and 
generation of the cosmological baryon asymmetry. We have shown conclusively (up to mass 



effects [^, |T9[) that in a i? — L = universe (with B — L conserved) sphaleron erasure of 
a primordial asymmetry does indeed occur in the standard model. The final asymmetry in 
this case is exponentially sensitive to model parameters and it may be possible to preserve 
the asymmetry in this case in an extension of the standard model. 

Even in models which do not conserve B — L, sphaleron erasure of a baryon asymmetry 
may occur if other i?/L- violating interactions are included in any extension of the standard 
model. A simple example of such an interaction is the AL = 2 operator induced from neu- 
trino mass generation via the seesaw mechanism. The preservation of the baryon asymmetry 
places restrictions on such operators and thus on the neutrino mass matrix. In the literature 
several widely differing bounds on neutrino masses have been given. Here we have clarified 
these bounds and shown again the importance of the e^ decoupling. In particular we have 
shown that erasure may occur only between T^. and the electroweak phase transition tem- 
perature, Tc. We have also shown how the role of e^ decoupling should be applied to bounds 
on other B/L violating operators. 
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Finally, we have demonstrated that because of e^j decoupling, an initial lepton or e^j 
asymmetry can be used to generate a baryon asymmetry even when B = initially. In such 
cases, not only do we drastically weaken upper bounds on new operators, but even imply a 
lower bound such that they would be able to induce a nonzero B — L and ultimately the 
baryon asymmetry. 
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Figure Captions 

Fig. 1 (a) Higgs boson decay into riglit-lianded electron and lepton doublet; Higgs and 
gauge boson scattering of (b) L^ into Cr and (c) bosons into L^ and Cr, lepton number 
violating scatterings of (d) e^^" or e^z/^ into Higgs bosons and (e) e^-Higgs into (£ or 
z/£)-Higgs. 

Fig. 2 mH{T)/T as a function of the zero temperature Higgs boson mass mo for top quark 
masses of 180 GeV (top curve), 135 GeV (middle curve) and 90 GeV (lower curve). 

Fig. 3 (a) s-channel scattering of third-generation quark doublet and right-handed top an- 
tiquark into cr and L^; (b) t-channel scattering processes. 

Fig. 4 Diagrams contributing to cr production by Higgs and gauge boson scattering. The 
first two are present for both B and W gauge bosons, whereas the last exists only for 
B. 

Fig. 5 a) Contours of constant B final /Beq, when the scattering contribution to the chirality 
changing rate is omitted. Solid lines correspond to the present work; for comparison 
the previous result of CDEO is shown by the dashed lines, b) The solid lines are 
contours of constant Bjinai/Beq with the full interaction rate, including scatterings. 
Dashed lines are contours of constant thermal Higgs mass parameter xh- 

Fig. 6 a) Evolution of the asymmetries Le^, Lg^ and i? as a function of temperature, when 
electron lepton number is conserved at low temperatures, using the initial value i^e^.p = 
10~^. b) Evolution of absolute values of asymmetries with electron lepton number 
violation corresponding to {rriei) = 5 keV. B is initially negative and the dip in the 
corresponding curve marks the place where B becomes positive. Similarly Lg^ changes 
sign from positive to negative. Lg^ remains negative over the whole range. 
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Fig. 7 The shaded area between the curves C and D shows the region consistent with the 
observed baryon asymmetry for some choice of parameters rrit and itlh- Curve C 
corresponds to B,mt,mH = 4 x 10~^^, 90 GeV, 60 GeV and curve D to B^mt^niH = 
6 X 10-i\ 90 GeV, 1 TeV. 
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